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ABSTRACT 

A  physical  description  of  the  destabilizing 
mechanism  responsible  for  certain  pilot  induced 
lateral-directional  instabilities  is  included 
herein.   In  addition,  the  theoretical  and  experi- 
mental effort  propcrting  to  support  the  physical 
arguments  is  included. 

The  basis  for  this  study  was  a  pilot  induced 
instability  known  to  exist  for  the  X-15  research 
aircraft-pilot  combination  for  certain  conditions 
of  flight.   Accordingly,  a  major  portion  of  the 
effort  was  concentrated  upon  the  analysis  of  a 
specific  known  instability.   Nevertheless,  it  is 
the  opinion  of  the  authors  that  the  resultant 
conclusions  are  quite  general. 


SECTION  I 
INTRODUCTION 

The  investigation  upon  which  this  report  is  "based 
was  concerned  with  the  study  of  a  type  of  pilot  induced 
instability  of  the  lateral-directional  short  period  mode. 

In  the  more  usual  case,  an  instability  of  this 
nature  can  be  explained  as  a  function  of  the  character- 
istic response  time  of  a  human  being,  in  which  case  the 
problem  is  quite  simply  that  the  human  controller  can- 
not react  to  the  error  signals  with  sufficient  speeds 
Quite  obviously,  a  problem  of  this  nature  is  closely 
related  to  the  frequency  of  the  mode  of  motion. 

Some  recent  studies,  however,  have  indicated  the 
possibility  that  certain  pilot  induced  -  or  closed 
loop  -  instabilities  can  occur  which  do  not  appear  to 
be  strictly  a  function  of  frequency.   A  study  by 
Askenas  and  McRuer,  (Ref.  1)  based  upon  root  locus 
analyses  of  a  pilot-airplane  combination,  showed  the 
possibility  of  an  instability  of  this  type.   This 
study  further  indicated  certain  combinations  of  the 
aerodynamic  stability  derivatives  believed  to  be 
unfavorable. 

More  recently,  an  analysis  was  made  by  Taylor 
(Ref.  2)  as  a  result  of  simulator  studies  of  the  X-15 
research  aircraft.   These  simulator  studies  predicted 
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certain  flight  condition  boundaries  within  which  the 
X-15  closed  loop  system  would  exhibit  an  unstable 
lateral  directional  short  period  mode.   (Dutch  roll 
mode).   That  this  region  did  in  fact  exist  was  verified 
by  flight  tests. 

Both  of  these  studies  relied  heavily  on  the  root 
locus  technique  for  an  explanation  of  how  the  magnitude 
of  pilot  initiated  control  movements  could  destabilize 
an  otherwise  stable  system.   Consequently,  although  the 
studies  showed  that  the  pilot  could  destabilize  the 
system,  the  actual  destabilizing  mechanism  was  not 
explained.   This  study  was  initiated  in  an  effort  to 
provide  this  explanation. 

This  study  was  conducted  in  three  separate,  yet 
complementary,  phases.   The  theoretical  phase  consisted 
of  an  examination  of  the  basic  equations  of  motion  and 
the  significance  of  certain  parameters  upon  the  resultant 
dynamic  behavior  or  a  given  pilot-airplane  system.   Root 
locus  analyses  formed  a  major  portion  of  this  phase. 

What  has  been  called  the  computer-simulator  phase 
was  essentially  a  study  of  the  same  parameters  investi- 
gated theoretically,  but  with  the  aid  of  an  electronic 
analog  computer.   In  addition,  the  analog  computer  was 
used  as  a  fixed  base  simulator  for  the  investigation  of 
an  actual  closed-loop  system. 
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The  flight  phase  consisted  of  flight  tests  in- 
corporating the  variable  stability  N Avion  aircraft 
described  herein.   This  phase  was  primarily  initiated 
in  an  effort  to  provide  as  realistic  a  simulation  as 
possible  of  an  unstable  pilot-airplane  combination. 

In  general,  it  may  be  stated  that  the  results  of 
all  three  phases  contributed  significantly  to  the  con- 
clusions which  resulted  from  this  investigation. 
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SBCTIOK  II 
EQUIP  --JL:T 

The  variable  stability  flight  simulator  used  in 
this  investigation  was  a  North  American  NAvion  airplane. 
To  provide  the  variable  stability  capability,  a  modified 
Minneapolis-Honeywell  autopilot,  USA.F  Type  E-12,  had 
been  installed.   Modifications  to  the  autopilot  consist- 
ed of  the  addition  of  supplementary  rate  gyros  and  a 
sideslip  transducer.   Signals  from  the  transducers t 
proportional  to  the  measured  quantities,  were  introduced 
as  input  signals  to  the  autopilot,  resulting  in  control 
deflections  proportional  to  the  measured  quantity.   This 
made  the  system  capable  of  modifying  the  airplane  stabili- 
ty derivatives.   Data  collection  and  recording  were 
accomplished  by  means  of  an  A3C0P  pulse  width,  frequency 
modulated  telemetering  system.   An  analog  computer  was 
used  in  the  theoretical  development  of  the  problem,  and, 
in  conjunction  with  a  control  stick-operated  potenti- 
ometer, provided  fixed-base  simulation. 

TEST  AIRPLANE 

The  NAvion  is  an  all  metal,  low  wing,  single-engined 
monoplane.   The  engine,  a  Continental  E-185 ,  drives  a 
variable  pitch  Hartzel  propeller  and  is  rated  at  185 
horsepower  -  maximum  continuous  at  sea  level  at  2300  rpm. 
The  control  surfaces  are  of  conventional  design.   A 
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strearalined  static  balance  is  fixed  at  the  outboard  end 
of  each  aileron.   The  trim  tabs  for  the  aileron  and 
rudder  are  of  the  fixed-bend  type,  and  the  elevator 
trim  tab  is  adjustable  from  the  cockpit.   The  physical 
characteristics  of  the  NAvion  are  listed  in  Table  I. 
A  photograph  of  the  test  airplane  is  shown  in  Fig.  1. 
The  airplane  was  equipped  for  this  investigation  with  a 
three-axis  variable  feedback  autopilot. 

AUTOPILOT 

The  basic  autopilot  is  capable  of  altering  the 
stability  of  the  airplane  in  the  pitch,  roll  and  yaw 
axes  through  three  channels  of  operation.   For  this 
investigation,  the  pitch  channel  was  rendered  inoper- 
ative, with  the  exception  of  an  unmodified  elevator 
control.   Bach  channel  is  composed  of  an  AC  series 
summing  network  which  combines  signal  voltages  from  the 
sensing  transducers,  the  electrical  flight  controls  and 
the  knob  actuated  trim  controller.   The  difference  be- 
tween this  combined  signal  voltage  and  a  voltage  feed- 
back signal  from  the  servo-drum  position  transducer,  an 
error  signal,  is  applied  to  an  amplifier,  phase  dis- 
criminated and  transferred  into  servo  drum  rotation  by 
means  of  dual  power  relays.   The  servo  drum  rotation 
moves  the  control  system  of  the  aircraft.   The  motion  of 
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the  control  surface  continues  until  error  signal  voltage 
at  the  input  to  the  amplifier  is  reduced  to  a  value 
below  the  threshold  of  the  system.   This  value  is  about 
fifty  millivolts. 

The  basic  autopilot  had  been  modified  by  deleting 
some  of  the  original  provisions  and  by  supplying  ad- 
ditional feedback  loops.   The  directional  coupler  feed- 
back loop,  the  flight  control  stick  trim  provisions,  and 
the  coordinated  turn  provisions  had  been  eliminated. 
Sideslip. angle  and  roll  rate  feedback  loops  and  a  cross- 
control  signal  proportional  to  aileron  control  stick 
deflection  had  been  added  to  the  yaw  channel.   Sideslip 
angle  and  yaw  rate  feedback  loops  had  been  added  to  the 
roll  channel.   The  manually  p.djustable  ratio  potentiometers 
were  relocated  on  a  gain  control  panel  placed  between 
the  pilots*  seats.   An  electronic  safety  device  was 
installed  to  prevent  "hard-over"  signals  by  automatically 

disengaging  the  autopilot  if  an  error  voltage  exceeding 

r 

approximately  3  volts  should  instantaneously  appear  at 
the  amplifier  input. 

Quantities  to  be  telemetered  were  taken  from  the 
appropriate  terminals  on  the  main  terminal  board  to  a  15 
channel  filter-conversion  unit  which  converted  the  voltage 
to  DC.   A  gain  control  for  each  channel  permitted  the 
DC  voltage  to  be  adjusted  to  that  required  by  the  tele- 
metering system  (0-5  volts).   Prom  the  filter-conversion 
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unit,  the  DC  voltages  representing  the  measured  quanti- 
ties were  taken  to  an  auxiliary  terminal  strip  and 
connected  to  the  telemeter  transmitter  unit. 

The  electrical  controls  are  located  in  the  right 
seat  of  the  NAvion,  and  the  standard  manual  controls  in 
the  left  seat.   Reference  3  described  a  spring  centering 
device  and  fluid  dashpots  which  had  been  incorporated 
in  the  electric  stick.   These  had  provided  objectionable 
feel  characteristics. and,  according  to  the  recommendations 
in  Reference  3,  were  removed  prior  to  the  flight  test 
phase  of  this  investigation. 

INSTRUCT  AT  I  OH 

The  transducers  in  the  feedback  loops  of  the  modified 
autopilot  provided  signals  to  the  autopilot  amplifier 
through  the  AC  summing  networks  and  provided  the  signals 
required  to  measure  the  aircraft  motions.   The  flight 
condition  data  (airspeed,  altitude,  and  outside  air 
temperature)  were  obtained  from  standard  aircraft 
instruments.   Errors  in  these  instruments  were  assumed 
negligible.   The  physical  characteristics  of  the  com- 
ponents of  the  instrumentation  system  are  listed  in 
Table  II. 

Two  of  the  modified  rate  gyros  were  located  on  the 
autopilot  chassis  on  the  main  equipment  table  and  two 
were  mounted  on  the  floor  of  the  aircraft  in  the  equip- 
ment compartment o 
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The  use  of  a  single  sideslip  vane  was  possible 
because  it  contained  two  separate  pot intiometers.   The 
sideslip  vane  was  installed  on  a  boom  extending  four 
feet  ahead  of  the  wing  leading  edge  at  the  wing  tipv 
It  was  assumed  that  this  distance  was  adequate  to  mini- 
mize measurement  errors  due  to  the  wing  pressure  field. 
The  position  of  the  sideslip  vane  was  such  that  the  yaw- 
ing rate  of  the  aircraft  about  the  vertical  body  axis 
would  affect  the  sideslip  angle  measurement  slightly, 
but  this  effect  was  neglected. 

The  telemeter  transmitter  unit  consists  of  a 
rotary  switch  sampling  each  of  43  input  channels  plus 
two  synchronizing  channels  at  the  rate  of  20  times  per 
second.   The  sampled  data  are  converted  to  pulse  width 
form  by  a  keyer  unit  and  transmitted  as  a  UHF  frequency 
modulated  signal  to  the  telemeter  ground  station.   The 
filter-conversion  unit  limited  the  data  measurements  to 
14  quantities  plus  a  full  scale  reference  voltage.   By 
jumper  wiring  on  the  auxiliary  telemetering  terminal 
strip  it  was  possible  to  sample  a  given  quantity  more 
than  once  for  each  revolution  of  the  switch.   Specif i- 

i 

cations  of  the  telemeter  transmitter  unit  are  listed  in 
Table  II. 

GROUIiD  STATION 

The  ground  station  data  collection  equipment  con- 
sisted of  an  A3C0P  li  series  PW/FM  ground  station,  an 
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Ampex  Model  309C  dual  track  tape  recorder  and  associated 
graphical  recorders. 

The  ASCOP  ground  station  received  the  modulated 
signals  from  the  airborne  unit,  demodulated  and  decoded 
the  signals  to  provide  for  each  of  the  43  channels  of 
information,  and  provided  a  continuous  voltage  output 
to  the  recorders  representing  the  measured  flight  data. 

The  tape  recorder  had  provisions  for  simultaneous 
recording  of  telemeter  and  voice  transmissions.   This 
permitted  the  recording  of  description  of  runs,  pilot's 
comments  and  ground  station  operator's  comments  along 
with  the  telemeter  data. 

Sanborn  four  channel  pen  recorders,  Model  154-lOOB, 
were  used  in  the  graphical  recording  phases. 

ANALOG  COMPUTER 

The  analog  computer  used  in  this  investigation  was 
a  Goodyear  Aircraft  Corporation  Model  L3  (G-EDA)  linear 
electronic  differential  analyzer.   Twenty-four  auto- 
matically stabilized  DC  computing  amplifiers  were  avail- 

7 
able  with  open-loop  gain  greater  than  5  x  10  ,  and  of 

negligible  drift.   The  computer  incorporated  an  auto- 
matic error  indicator  and  had  a  guaranteed  accuracy  of 
one  percent.   Provisions  were  available  for  accurately 
setting  computer  board  potentiometers  by  the  use  of  a 
special  calibration  potentiometer  and  null  indicator. 
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At tempts  to  "fly"  the  computer  were  conducted  by 
means  of  the  simulator  shown  in  Figure  2.   This  consisted 
of  a  conventional  control  stick  which,  when  deflected, 
varied  the  voltage  at  the  <^a  terminals  of  the  analog 
computer  circuit  by  changing  the  setting  of  a  potenti- 
ometer. 
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SECTIOW  III 

PROCEDURE  r 

A,   In-Flight  Simulation 

In-flight  simulation  of  the  pilot  induced  lateral- 
directional  oscillation  was  accomplished  with  the 
variable  stability  TTAvion.   It  was  expected  that  results 
obtained  in  this  manner  would  provide  a  more  suitable 
basis  for  a  physical  explanation  of  the  phenomenon  than 
results  obtained  from  a  fixed  base  simulator,  since  the 
full  range  of  visual  and  kinesthetic  motion  cues  would 
be  present. 

Reference  2  points  out  uncontrollable  combinations 
of  Mach  number  and  angle  of  attack  encountered  with  the 
X-15  research  aircraft  with  lateral  stability  augment- 
ation off  and  the  pilot  attempting  to  control  bank  angle 
in  a  normal  manner  with  ailerons  alone.   Using  these 
conditions  as  a  guide,  values  were  obtained  for  the 
dimensional  stability  derivatives  to  be  reproduced  in 
the  K Avion.   These  are  listed  in  Table  III.   The  calcu- 
lation procedure  is  shown  in  Appendix  A.   These  deriva- 

f 

tives  were  dimensionalized  for  a  dynamic  pressure  of 
153  lbs. /ft.  ,  corresponding  to  one  g  flight.   This 
choice  of  a  relatively  low  dynamic  pressure  was  made  in 
order  that  the  capabilities  of  the  autopilot/airplane 
system  not  be  exceeded;  it  was  anticipated  that  the 
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frequency  of  the  Dutch  roll  oscillations  would  be 
reduced,  while  the  character  of  the  pilot-induced 
instability  would  be  retained.   Simulation  at  higher 
dynamic  pressure  was  confined  to  the  fixed  base  simulator. 

Dimensional  derivatives  for  the  NAvion  were 
obtained  for  6500  ft.  standard  density  altitude  and 
120  miles  per  hour  true  airspeed  (Ref.  4).   Autopilot 
feedback  gain  constants  were  then  determined  which  would 
alter  these  derivatives,  making  them  equal  to  those 
chosen  as  typical  for  the  X-15.   This  procedure  is  shown 
in  Appendix  A. 

feedback  gain  constants  of  the  autopilot  could  be 
altered  by  adjustment  of  potentiometer  settings  in  the 
cockpit  of  the  NAvion.   In  order  to  determine  gain 
potentiometer  settings  which  would  provide  the  required 
constants  it  v/as  necessary  to  calibrate  the  system.' 

Each  autopilot  rate  gyro  was  calibrated  by  removing 
it  from  the  airplane  and  remounting  it  on  an  electrically 
driven  turntable  which  could  be  rotated  at  a  controllable 
rate.   In  this  manner  control  surface  deflections  pro- 
duced by  given  angular  velocities  could  be  established. 
The  deflections  were  measured  by  a  propeller  protractor 
for  various  gain  potentiometer  settings. 

The  sideslip  vane  was  calibrated  by  moving  it  through 
a  specified  angle  and  simultaneously  measuring  the  control 
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surface  deflections  produced  at  different  gain  potenti- 
ometer settings.   Electric  stick  calibration  was  obtain- 
ed in  the  same  manner,  a  propeller  protractor  being 
used  to  measure  stick  deflection. 

The  results  of  the  calibration  are  shown  in  Figure 
3,  as  curves  of  Feedback  Gain  Constant  versus  Gain 
Potentiometer  Faceplate  Setting.   The  numbering  of  the 
constants  corresponds  to  the  number  of  the  potentiometer 
by  means  of  which  each  may  be  adjusted.   The  cockpit 
panel  arrangement  is  shown  in  Figure  4. 

In  order  that  a  quantitative  interpretation  of  the 
telemetered  motion  histories  of  the  aircraft  could  be 
made,  it  was  also  necessary  to  calibrate  the  system 
for  telemetry.   This  was  performed  in  exactly  the  same 
manner  as  the  calibration  described  above,  with  the 
exception  that  the  measured  output  quantity  in  this 
case  was  voltage  at  the  telemetry  transmitter  package, 
rather  than  control  surface  deflection.   It  was  decided 
that  the  following  quantities  would  provide  sufficient 
information  for  the  analysis  to  be  performed:   roll  rate, 
yaw  rate,  stick  deflection  and  sideslip  angle.   The 
results  of  the  telemetry  calibration  are  shown  in 
Figure  5  as  curves  of  tne  foregoing  quantities  versus 
Telemetry  Input  Voltage. 
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Longitudinal  feedbacks,  (pitch  angle,  pitch  rate, 
velocity  and  angle  of  attack)  were  not  employed  in  this 
investigation.   The  elevator  gain  potentiometer  was 
adjusted  to  a  value  which  gave  a  satisfactory  elevator 
response. 

Certain  deficiencies  existed  in  the  equipment  which 
made  it  difficult  to  achieve  an  exact  simulation.   Among 
these  were  control  cable  stretch,  structural  elasticity 
of  the  aircraft,  inaccurate  resolution  due  to  finite 
autopilot  "peck  sizes,"  very  small  potentiometsr  settings 
and  the  difficulty  of  accurately  reproducing  a  potenti- 
ometer setting.  As  a  result  of  these  deficiencies  it 
was  found  that  the  simulation  procedure  outlined  in  this 
section  provided  only  a  first  approximation  to  the  desired 
aircraft  response.   It  was  necessary,  therefore,  to 
further  adjust  the  cockpit  gain  potentiometer  settings 
in  order  to  achieve  the  desired  results. 

This  arbitrary  correction  to  the  calculated  potenti- 
ometer settings  would  necessitate  matching  of  aircraft 
responses  to  analog  computer  responses  if  the  objective 
was  exact  simulation.   A  reasonably  good  qualitative 
simulation  of  such  motions  can  be  readily  achieved, 
however,  and  was  considered  to  be  sufficient  for  the 
purpose  of  this  investigation. 
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B.   Analog  Computer  Study 

Three  basic  tools  were  used  to  accomplish  the 
computer  phase  of  this  investigation.   These  were: 

1.  the  (PTDA  Analog  Computer t 

2.  the  Sanborn  Recorder,  and 

3.  the  Control  Stick  Simulator. 

The  computer  and  recorder  were  used  to  do  the  open  loop 
studies  of  the  system,  while  all  three  units  were  used 
to  complete  the  closed  loop  work. 

The  equations  mechanized  in  the  computer  are  listed 
and  discussed  in  detail  later  in  this  report. 

When  open  loop  response  to  step  inputs  by  the 
computer  matched  calculated  values  as  to  time  to  damp 
to  1/2  amplitude  and  period  of  oscillation,  attempts 
were  made  to  "fly"  the  computer.   These  closed  loop 
trials  were  accomplished  by  using  the  control  stick 
simulator  in  the  same  manner  as  the  control  stick  in  an 
aircraft.   The  time  histories  of  bank  angle,  roll  rate, 
sideslip  angle,  and  yaw  rate  on  the  recorder  were  used 
as  references.   The  control  stick  simulator  produced 
input  voltages  to  the  computer  analogous  to  aileron 
inputs  to  the  basic  airplane. 

Three  separate  configurations  of  the  basic  equations 
of  motion  were  investigated.   These  were: 
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A.  O^o =  0  and  low  dynamic  pressure  (153  psf). 

B.  C*0  =  10°  and  low  dynamic  pressure  (153  psf), 

C.  0Co  =  10°  and  high  dynamic  pressure  (1000  psf). 

The  wiring  diagrams  for  the  computer  study  are 
shown  in  Figures  6  and  7. 
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section  iv 

discussion  ai;d  results 

Theoretical  Development; 

During  the  past  several  years,  considerable  emphasis 
has  been  placed  on  the  determination  of  what  parameters 
can  best  be  used  to  measure  the  lateral-directional 
handling  qualities  of  an  airplane.   As  a  result,  a  consider- 
able body  of  literature  has  come  into  being  and  certain 
airplane  response  characteristics  have  come  to  be  re- 
garded as  unfavorable  or  undesirable  from  the  standpoint 
of  control  by  a  human  pilot.   (Ref.  1).   Primarily  as  a 
result  of  the  human  pilot  being  involved  as  a  part  of  the 
feedback  control  loop,  investigations  of  this  nature  have 
not  been  easy  to  interpret  either  physically  or  mathe- 
matically. 

This  particular  investigation  has  been  concerned 
with  a  rather  special  segment  of  the  lateral-directional 
handling  qualities  problem  -  specifically  the  case  where 
the  inclusion  of  a  pilot  in  the  control  loop  has  resulted 
in  an  unstable  Dutch  roll  oscillation,  although  the  open 
loop  Dutch  roll  response  of  the  aircraft  was  known  to  be 
stable.   In  short,  this  investigation  has  been  concerned 
with  the  specific  problem  of  a  pilot  induced  Dutch  roll 
instability. 

The  problem  of  a  pilot-airplane  combination  of  the 
nature  described  above  could  have  been  regarded  as  general 
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and  an  overall  solution  sought.   However,  it  was  felt 
by  the  authors  that  a  more  concentrated  effort  could  be 
made  if  the  problem  were  approached  from  the  opposite 
point  of  view  -  namely  a  detailed  theoretical  and 
experimental  investigation  using,  as  a  guide,  a  pilot- 
airplane  combination  for  which  the  results  were  known. 
The  reasoning  followed  was  that  a  detailed  study  of  a 
particular  situation  might  then  yield  a  theory  which 
could  be  investigated  in  a  more  general  sense. 

According  to  Ref.  2,  the  X-15  research  airplane- 
pilot  combination  was  found  to  result  in  an  unstable 
lateral-directional  oscillation  for  certain  flight 
conditions  if  the  stability  augmentation  system  was  made 
inoperative.   For  this  reason,  the  X-15  airplane-pilot 
combination  was  the  closed  loop  system  chosen  for  the 
detailed  study.   Specifically,  the  instability  associated 
with  a  flight  Maoh  number  of  3  and  a  ten  degree  angle  of 
attack  was  investigated. 

The  fact  that  the  above  configuration  v/as  the  one 
chosen  for  study  placed  several  important  restrictions 
upon  the  problem  and  thereby  upon  the  resultant  conclusions. 
First,  although  the  above  configuration  was  handled 
throughout  this  study  as  though  the  aircraft  had  been  in 
steady  flight,  in  reality  the  X-15  encountered  the  above 
stated  flight  condition  only  in  transitory  flight. 
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Accordingly,  the  aerodynamic  stability  derivatives  were 
continuously  changing,  with  the  result  being  that  the 
pilots  of  the  X-15  were  advised  to  refrain  from  using 
the  rudder  as  a  control.   Therefore,  this  entire  investi- 
gation was  conducted  assuming  the  ailerons  to  be  the  only 
type  of  lateral -directional  control  available  ,  although 
the  effect  of  using  rudder  is  discussed  in  the  flight 
phase  section  of  this  discussion.   A  second  condition 
implied  by  the  problem  as  stated  is  that  the  frequency  of 
the  oscillation  was  not  a  factor.   That  this  must  be  so 
can  be  reasoned  by  noting  that  the  frequency  is  more  a 
function  of  the  dynamic  pressure  than  of  Mach  number, 
while  the  instability  has  been  described  as  being  a 
function  of  Mach  number  and  angle  of  attack.   Thus,  the 
solution  sought  was  not  one  dependent  upon  the  inability 
of  a  human  being  to  respond  with  sufficient  speed  to  an 
error  signal.   In  effect,  a  suitable  explanation  to  the 
problem  encountered  by  the  X-15  would  have  to  be  valid 
even  if  the  pilot  were  assumed  to  respond  without  a  time 
lag. 

The  stability  derivatives  required  for  a  detailed 
study  of  the  above  condition  were  obtained  from  Ref .  5 
and  are  listed  in  Table  TV.   These  derivatives  lead  to 
autopilot  potentiometer  gain  constants  listed  in  Table  V. 
It  is  noteworthy  that  only  two  of  the  stability  derivatives 


-28- 
were  significantly  different  from  what  one  might  ordi- 
narily expect.   Namely,  bothCn(£q,  (aileron  yaw)  and  Cfla 
(dihedral  effect)  had  a  positive  sign.   Thus,  these  two 
derivatives  immediately  became  suspect  and  it  was  expected 
that  a  theoretical  analysis  of  the  problem  would  show 
either  one  or  both  of  these  derivatives  to  be  involved 
in  an  explanation  of  the  instability.   In  this  rather 
special  case  of  a  single  control,  pilot-induced  insta- 
bility, the  unusual  sign  of  the  aileron  yaw  became  of 
particular  interest. 

A  proposed  theoretical  explanation  to  the  problem 
encountered  by  the  X-15  aircraft  was  published  by  Taylor 
(Ref.  2).   In  his  work,  Taylor  showed,  by  root  locus 
techniques,  that  the  Dutch  roll  mode  of  an  aircraft  with 
the  characteristics  exhibited  by  the  X-15  could  be  made 
unstable  if  a  pilot  gain  proportional  to  roll  angle  and 
roll  rate  was  included  in  the  feedback  loop.  According 
to  Taylor,  the  ratio  of  ^/a)^  could  be  used  to  define 
the  region  where  an  instability  might  occur  -  a  ratio 
greater  than  one  being  indicative  of  an  unfavorable 
condition.   The  human  transfer  function  used  by  Taylor 
in  his  analysis  was  determined  empirically  and  found  to 
be  K(l  +  .57s).   Thus  a  pilot  response  proportional  to 
both  roll  angle  and  roll  rate  was  assumed.   The  analysis 
followed  by  Taylor  was  quite  believable  for  two  reasons. 
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First,  the  transfer  function  was  seen  to  be  independent  of 
a  pilot  time  lag  which  would  seem  to  isolate  any  effect  of 
frequency  from  the  basic  control  problem.   Secondly,  the 
analysis  as  carried  forth  by  Taylor  resulted  in  the  con- 
clusion that  the  basic  cause  of  the  problem  was  the  positive 
sense  of  the  dihedral  effect  (CjJa).   What  Taylor's  analysis 
did  not  do  was  present  a  physical  explanation  of  the  phe- 
nomenon. 

In  an  effort  to  assure  that  the  problem  had  been 
interpreted  correctly  by  the  authors,  the  initial  theo- 
retical effort  was  directed  toward  matching  the  results 
of  Taylor's  investigation.   Thus  the  human  transfer 
function  proposed  by  Taylor  was  used  and  a  root  locus 
study  of  the  effect  of  varying  pilot  gain  undertaken. 
The  set  of  equations  used  was  as  follows: 
-£-  d^  J>  &a 

(1)   yP 

o 

It  is  worthy  of  note  that  the  above  set  is  for  a  prin- 
cipal axes  coordinate  system.   This  particular  coordinate 
system  was  chosen  for  this  study  inasmuch  as  the  stability 
derivatives  as  obtained  from  Ref .  5  had  been  calculated 
with  respect  to  this  system  of  axes. 

The  result  of  the  root  locus  study  can  best  be 
explained  by  considering  separately  the  effect  of  varying 
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the  pilot  gain,  K,  and  the  trim  angle  of  attack,  (V0. 
The  fact  that  <X0  was  a  significant  parameter  actually 
was  determined  during  the  analog  computer  phase  of  this 
program  and  is  discussed  in  detail  in  connection  with 
that  part  of  the  program.   Here  the  fact  that  Oi0 is  import- 
ant will  merely  be  stated  and  then  substantiated  by  the 
root  locus  technique. 

Assuming  the  pilot  gain  to  be  the  variable  of  interest, 
the  determinant  of  the  coefficients  of  (1)  was  expanded 
with  the  following  root  locus  equation  resulting: 

(2)     (^^rXA-XsH  AZ  +   atyCOnyA+COny)    uClW     n 
(\z-   ZtfiUnfiK  +  U>np*X\  +   -57  A)  HX~ 

The  derivation  of  the  above  equation  is  included  in 

Appendix  B. 

Figure  8a  was  the  type  of  root  locus  which  resulted 

from  the  above  equation  when  OCo  was  assumed  to  be  negligible*. 
This  root  locus  showed  the  inclusion  of  a  proportional 
pilot  gain  to  be  a  stabilizing  factor  at  zero  angle  of 
attack,  provided  the  stability  derivatives  could  be 
assumed  constant.   Furthermore,  a  study  of  the  parameters 
which  determine  the  value  of  C0n£>  (^q.  5  in  Appendix  B), 
reveals  that  an  increase  in  either  of  the  two  suspect 
stability  derivatives  (Cn£  and  Qflft)  would  actually  de- 
crease the  value  of Cd^  and  thus  tend  to  stabilize  the 
closed  loop  system.   This  is  a  result  of  the  signs  of 
both  derivatives  being  reversed  from  what  might  be  termed 
normal. 


imag .  axis 


real  axis 
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Fig.    8a.    Roo+  Locus  for  cx>0  ,  0<K<°° 


The  next  step  was  to  include  the  correct  value  of  0£0 
(10  degrees).   Prom  equations  (1)  it  can  be  seen  that  the 
inclusion  of  0t0  has  no  effect  upon  the  denominator  term 
of  the  root  locus  equation  but  does  alter  the  character- 
istic equation  somewhat.   The  added  terms  which  are  the 
result  of  including  OC0  are: 


(3)  <Xf 


^(cy^-c^cOA-^M 


It  can  be  seen  from  equation  (3)  that  with  a  positive 
Cj^and  all  other  derivatives  of  the  usual  sign,  the 
effect  of  including  (X0is  to  decrease  the  A  and  A 
terms  of  the  characteristic  equation.   This  had  little 
effect  upon  the  roll  and  spiral  modes  but  did  result  In 
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both  the   damping  ratio   (§m  )   and   the  undamped  natural 
frequency   (00rn|> )   of  the  Dutch  roll  mode  being  appreci- 
ably decreased.      The  root   locus  which  resulted  is  depicted 
by  Pig.    8b. 


imacj  axis 


real   axis 


Fig  8b.     Roo+  Locus  for   OCa=  10°,     0<K<oO 


This  is  similar  to  what  Taylor  obtained.   Furthermore, 
it  can  be  seen  from  equation  (3)  that  the  destabilizing 
effect  of  0CQ   can  be  attributed  to  the  positive  sign  of 
the  dihedral  effect.  (Cfia). 

Having  shown  the  importance  of  the  trim  angle  of 
attack,  it  was  deduced  that  a  root  locus  using  trim 
angle  of  attack  (  (X0 )  as  the  varying  quantity  might  be 
informative.   This  study  was  carried  out  by  assuming  the 
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stability  derivatives  to  "be  independent  of  the  trim 
angle  of  attack.   Thus  only  the  orientation  of  the  air- 
craft with  respect  to  its  flight  path  was  assumed  to  be 
variable.   The  root  locus  equation  for  this  investigation 
had  the  following  form: 

(4)(XAr)(X+Xs)(X2"+2^COn<P->-COn(pZ)  .  k,      - 

A(X  +  K'/K')         ^a& 

with  the  denominator  terms  found  from  equation  (3)  and  with 

Figure  8c  shows  the  root  locus  for  0^°^o^00f  or  no 
pilot  gain  while  figure  8d  shows  the  same  result  for  an 
arbitrary  pilot  gain.   Calculations  for  the  conditions 
investigated  showed  the  system  with  pilot  gain  to  go 
unstable  for  a  trim  angle  of  attack  on  the  order  of 
twenty  five  percent  of  what  was  required  for  instability 
without  the  pilot  in  the  loop.   Conversely  it  was  noted 
that  for  a  trim  angle  of  attack  of  ten  degrees,  the  open 
loop  system  was  so  slightly  stable  that  a  small  pilot 
gain  could  be  expected  to  be  sufficient  to  cause  the 
closed  loop  instability.   Prom  equation  (4)  it  can  be 
seen  that  a  negative  Cjta  would  change  the  sign  of  K* 
and  in  turn  change  the  angle  condition  for  the  root  locus 
diagram.   For  this  case,  the  effect  of  increasing  OC© 
would  not  be  destabilizing. 
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Several  interesting  o"bservations  resulted  from  a 
study  of  the  root  loci  "briefly  discussed  above.  First, 
it  was  observed  that  if  the  signs  of  both  the  dihedral 
effect  (  CjJ/3  )  and  the  aileron  yaw  (Cn&a  )  nad  been' 
reversed,  there  would  have  been  virtually  no  change  in 
the  closed  loop  response  at  low  angles  of  attack.  How- 
ever, a  study  of  the  cof actor  of  the  (CX0)  terra  showed 
that  a  beneficial  effect  would  have  been  realized  as 
the  trim  angle  was  increased.   Secondly,  the  observation 
was  made  that  for  a  positive  dihedral  effect,  the  pro- 
verse  sense  of  the  aileron  yaw  definitely  tended  to 
improve  the  closed  loop  stability  while  an  adverse 
aileron  yaw  would  have  had  the  opposite  effect.   That 
is,  an  adverse  aileron  yaw  would  have  placed  the  complex 
zero  (Fig.  8b)  even  further  above  the  complex  pole. 
This  was  believed  to  be  particularly  significant  inas- 
much as  with  pilot  time  lag  removed  as  a  significant 
contributor  to  the  instability,  one  might  reason  that 
the  only  way  that  the  pilot  could  be  contributory  as  a 
destabilizing  factor  would  be  as  a  means  for  introducing 
an  unfavorable  forcing  function  while  actually  responding 
correctly  to  the  error  signal.   In  the  case  of  a  single 
control,  this  meant  an  unfavorable  input  into  either  the 
rolling  or  yawing  moment  equations  by  means  of  an  aileron 
deflection.   For  the  particular  case  investigated,  the 
root  locus  plots  showed  a  proverse  aileron  yaw  (Cn&a  ) 
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to  be  favorable  while  an  increase  in  the  trim  angle  (CX0) 
was  seen  to  have  a  destabilizing  effect.   Thus  it  was 
reasoned  that  the  probable  physical  effect  was  that  an 
increase  in  trim  angle  of  attack  opposed  the  proverse 
aileron  yaw  effect  until,  at  some  angle  of  attack,  the 
sense  of  the  sideslip  resultant  from  an  ailsron  deflection 
was  probably  reversed.   This  might  be  explained  by  noting 
that  a  simple  uncoupled  rolling  motion  about  the  body 
axis  of  an  aircraft  will  result, as  the  roll  angle  increases, 
in  the  conversion  of  angle  of  attack  into  sideslip  angle. 
For  example,  a  positive  roll  would  tend  to  produce  a 
positive  sideslip.   This  is  the  opposite  of  the  effect 
produced  by  a  proverse  aileron  yaw.   Thus  it  is  reason- 
able to  assume  that  at  some  point  the  effect  of  angle  of 
attack  and  of  proverse  aileron  yaw  might  be  equal.   At 
this  point  and  beyond,  an  increase  in  the  pilot  gain 
would  cause  the  pilot  to  introduce  a  destabilizing  force. 
Looked  at  in  the  light  of  Taylor* s  analysis,  this  point 
where  the  aileron  yawing  moment  due  to  aileron  deflection 
reversed  direction  is  probably  coincident  with  the  point 
where  ^t/LOy  becomes  greater  than  unity. 

In  summary,  the  theoretical  effect  agreed  quite  well 
with  Taylor's  analysis  although  it  is  felt  by  the  authors 
that  the  approach  involving  the  variation  of  trim  angle 
of  attack  provided  a  basis  for  the  physical  reasoning 
explained  above.   However,  it  can  be  seen  that  even  if 
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the  physical  effect  of  increasing  the  trim  angle  is  as 
stated,  the  mechanism  by  which  the  instability  occurs 
cannot  be  explained  from  the  theory  alone.   Fortunately, 
a  study  of  the  time  histories  of  both  the  analog  computer 
simulator  and  the  variable  stability  NAvion  did  provide 
an  insight  into  the  mechanism  involved  and  did  tend  to 
agree  with  the  hypothesis  suggested  by  the  theory. 

Analog  Computer  -  Simulator  Phase; 

As  mentioned  previously,  the  approach  used  in  this 
investigation  was  to  make  a  detailed  study  of  an  airplane- 
pilot  combination  which  results  had  shown  to  be  unstable. 
For  the  computer  phase  of  the  investigation  two  variations, 
along  with  the  basic  condition,  were  examined.   The  basic 
condition  was  the  X-15  research  airplane  flying  at 
M  =  3.0  with  an  angle  of  attack  of  ten  degrees.   The  two 
variations  were:   to  keep  all  other  conditions  the  same 
but  decrease  the  trim  angle  of • attack  to  zero,  and  to 
keep  all  other  conditions  the  same  while  increasing  the 
dynamic  pressure  from  153  to  1,000  pounds  per  square  foot. 
All  three  of  these  conditions  wiip.  be  discussed  below. 

Condition  I:   M  =  3.0,  q  =  153#/ft2.   (Xo  =  0 
For  this  condition,  the  basic  equations  of  motion, 
using  the  non-dimensional  X-15  coefficients,  reduce  to 
the  following  set  of  computer  equations: 
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100(3   =   -3.35(3   -    100^   +    1.0670 
.20    =      1.67(3+  .0345^ -,Q4£^  •+ 3.66<§q 
lOf    =  60(3-    .451  4*  -.04510  +  9.96a 

Both  step  aileron  inputs  and  "attempts-to-fly"  the 
computer  simulator  showed  this  condition  to  be  completely- 
manageable  without' the  slightest  tendency  to  go  unstable. 

This  fact,  in  conjunction  with  the  other  studies 
outlined  herein,  confirmed  the  significance  of  the  trim 
angle  of  attack  (  0CQ)   to  the  study  of  the  closed  loop 
instability. 

Condition  II;  M  =  5.0,  q  =  153#/ft2.  &o =  10  degrees 

This  condition  resulted  from  an  analysis  of  the 
results  of  Condition  I.   By  re-evaluation  of  the  basic 
equations,  using  0L0  =   10  degrees,  the  following  set  of 
equations  was  obtained: 

100/3  =  -3.35(3  -  1 00  4>  +-  1.067  0  +  17-45  0 
10  0  =  90(3  +  1 .725^  -  2.12  0  +-  83.9 So 
20<£    x         55(5    -    .902^  -.0902  0    +     9.08c§a 

These  equations  were  mechanized  into  the  computer  and 
resulted  in  the  predicted  pilot-induced  instability. 

A  step  aileron  input  and  an  "attempt-to-fly"  the 
computer  are  shown  in  Figures  9  and  10.   From  the  step 
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input  response  there  were  several  important  things  noted. 
First,  it  was  seen  that  the  motion  was  lightly  damped. 
Second,  there  was  considerable  excitation  of  the  Dutch 
roll  mode  from  a  coupling  of  the  rolling  and  yawing  modes 
of  motion.   Third,  although  the  sign  of  Cn£q  was  posi- 
tive -  leading  one  to  believe  the  yaw,  and  consequently 
the  sideslip,  due  to  aileron  should  have  "been  proverse  - 
the  total  effect  of  trie  parameters  affecting  the  side- 
slip was  just  the  opposite.   That  is,  the  sideslip  angle 
response  described  by  these  equations,  due  to  a  step 
aileron  input,  was  as  might  be  expected  from  an  adverse 
yawing  moment.   The  fourth  and  last  item  to  be  noted  was 
the  phase  angle  difference  between  the  angle  of  roll  and 
the  angle  of  sideslip.   It  could  be  seen  from  the  time 
histories,  that  the  two  motions  v/ere  almost  exactly  out 
of  phase.   The  significance  of  the  phasing  of  these  angles 
while  the  Dutch  roll  motion  was  in  progress,  as  well  as 
the  sign  of  the  sideslip  resulting  from  an  aileron  input, 
will  be  discussed  more  fully  later. 

Using  the  control  stick  simulator,  "attempts-to-fly" 
the  simulator  resulted  in  oscillations  of  increasing 
amplitude  if  the  "pilot"  responded  in  a  normal  manner 
to  bank  angle  and  roll  rate.   Using  the  stick  to  make 
aileron  inputs  in  an  effort  to  keep  the  wings  level 
resulted  in  an  increase  of  the  sideslip  ang]e  ,  due  to 
the  aileron  deflection.   Thus,  sideslip  grew,  fed  the 
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roll  excursions,  and  the  resultant  situation  was  one  of 

a  closed  loop  instability.   By  releasing  the  stick,  the 

oscillations  damped  out,  although  the  roll  angle  grew 

indefinitely  as  a  result  of  the  unstable  spiral  mode. 

An  "attempt-to-fly"  the  computer  by  using  the  /5 

technique  described  by  Petersen  in  Reference  6  was 

completely  successful.   The  oscillations  were  allowed  to 

build  up  using  normal  pilot  techniques,  and  then,  at 

some  point,  ailerons  were  used  primarily  to  stop  the 

r 
sideslipping,  and  secondarily  to  keep  the  wings  level. 

As  the  |3  trace  came  through  zero,  the  stick  controller 

would  give  the  stick  a  sharp  pulse  in  the  direction  the 

trace  was  moving.   In  the  airplane  this  would  correspond 

to  a  sharp  right  aileron  pulse  as  the  nose  of  the  airplane 

went  through  zero  degrees  sideslip  heading  toward  the 

right.   In  this  manner,  the  motion  could  be  damped  out 

within  a  few  cycles.   It  did  not  appear  feasible  to 

maintain  steady  state  conditions  using  this  technique, 

but  oscillations  resulting  from  normal  pilot  inputs  could 

be  rapidly  damped  out  in  this  manner.   Roll  angle  control 

was  maintained  by  varying  the  magnitude  of  the  pulse  and 

the  stick  position  betv/een  pulses. 

During  the  study  of  this  condition,  the  following 

observations  were  noted: 

1.  The  oscillatory  motion  was  lightly  damped, 

2.  There  was  considerable  coupling  between 
the  roll  and  sideslip  motions, 
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3.  The  sideslip  generated  by  an  aileron  deflection 
was  adverse  even  though Cn6a  was  positive. 

4.  During  Dutch  roll  oscillations , the  angles  of 
bank  and  sideslip  remained  approximately  180 
degrees  out  of  phase. 

5.  Normal  attempts-to-fly  were  unsuccessful,  and 

6.  The  j2>  technique  could  be  used  to  damp  out  any 
Dutch  roll  oscillations  that  developed. 

Condition  III:   M  =  3.0,  q  =  1000,f/ft2,  (Xo  =  10  degrees 

This  condition  was  investigated  primarily  to  see  if 
the  characteristics  of  the  motion  were  affected  by  changes 
in  dynamic  pressure.   The  basic  equations  of  motion  modi- 
fied to  increase  the  dynamic  pressure  were  written  as: 

20(3   =   -4.52(3    -  20  <£    +    1.440  +  3-490 
20    =    116.6(5 +2.32 <P    -2.840   +  l08.4<Sq 
49    =     71  .2(3-1. 2l2(f'-.  12120    +11. 72.  6Q 

These  equations  were  mechanized  in  the  computer  and 
both  step  aileron  inputs  and  nattempts-to-fly"  were 
observed.   The  basic  character  of  the  motion  did  not 
change.   All  the  comments  recorded  under  Condition  II 
are  valid  for  the  higher  dynamic  pressure.,  The  only 
significant  change  in  the  motion  was  the  frequency  of 
oscillation  of  the  Dutch  roll  mode.  Previously  the 
motion  had  a  period  of  about  five  seconds.   At  this  in- 
creased dynamic  pressure  the  period  was  reduced  to 
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approximritely  two  seconds.   This  did  not  alter  the  diffi- 
culty of  attempting  to  fly  the  configuration.   Nor  did 
the  increased  frequency  invalidate  the  fo  technique  of 
control. 

This  condition  can  be  summed  up  simply  "by  saying  it 
was  identical  to  Condition  II  except  things  happened  a 
little  faster.   However,  the  motion  was  not  so  fast  as 
to  make  the  controller^  task  too  difficult,  nor  so  fast 
that  the  character  of  the  motion  and  subsequent  responses 
to  the  motion  were  changed. 

One  result  of  the  above  observations  was  that  the 
frequency  of  the  motion  (within  the  limits  investigated) 
was  seen  to  be  an  insignificant  contributor  to  the  in- 
stability.  As  long  as  the  motion  proceeded  at  a  pace 
which  the  controller  could  follow,  and  as  long  as  the 
normal  control  motions  were  made,  the  basic  instability 
occurred.   Thus  the  frequency  of  the  motion  was  elimi- 
nated as  an  important  parameter  insofar  as  a  study  of 
the  basic  destabilizing  mechanism  was  concerned.   This 
confirmed  the  hypothesis  outlined  in  the  theoretical 
portion  of  the  discussion. 

Flight  Phase: 

The  flight  phase  of  this  investigation  involved  the 
use  of  the  variable  stability  KAvion.   The  physical 
characteristics  of  this  airplane,  its  simulation  equip- 
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ment  and  the  accompanying  instrumentation  are  described 
elsewhere  in  this  report. 

The  initial  flight  phase  was  "based  upon  an  in- 
stability encountered  by  the  X-15  research  aircraft  in 
a  steady  flight  (lift  =  weight)  condition  at  Mach  3  and 
at  ten  degrees  angle  of  attack.   The  dimensional  X-15 
derivatives  used  in  calculating  the  potentiometer 
settings  are  listed  in  Table  III.   A  direct  method  for 
simulating  the  X-15  side  force  equation  was  not  available, 
and  the  dynamics  introduced  by  the  autopilot  equipment 
were  not  known  by  the  authors,  hence  it  was  necessary  to 
correct  the  computed  autopilot  potentiometer  gain 
settings  in  flight  until  a  desired  condition  was  achieved. 

Flight  simulation  was  directed  toward  achieving  the 
desired  character  of  pilot  induced  instability,  rather 
than  toward  an  exact  simulation  of  the  dynamic  response 
of  the  X-15.  On   this  premise,  and  on  the  basis  of 
computer  studies,  it  was  decided  that  the  most  signifi- 
cant characteristics  to  be  simulated  were  the  following: 
1)   a  lightly  damped  Dutch  roll  mode  2)  a  roll  rate 
response  to  an  aileron  function  which  contained  a  signifi- 
cant Dutch  roll  effect. 

Cockpit  potentiometer  settings  were  adjusted  until 
the  two  above  characteristics  were  obtained.   It  was 
found  that  the  period  of  the  Dutch  roll  mode  obtained 
in  this  manner  was  approximately  three  seconds,  and  the 


-46- 

steady  state  roll  rate  of  the  aircraft  was  slightly  de- 
creased.  Thus  Dutch  roll  period  is  of  the  same  order 
of  magnitude  as  the  periods  of  the  Dutch  roll  modes 
studied  on  the  computer.   In  addition,  computer  studies 
indicated  that  the  frequency  of  the  Dutch  roll  mode  was 
unimportant  insofar  as  the  fundamental  problem  was 
concerned. 

Undoubtedly  the  most  important  characteristic  was 
simulated,  namely,  the  closed-loop  pilot-airplane  system 
exhibited  an  unstable  Dutch  roll  oscillation  in  spite  of 
the  fact  that  the  open  loop  response  was  lightly  damped. 
Thus,  as  in  the  case  of  the  X-15,  the  attempts  of  a 
pilot  to  improve  upon  a  marginal  condition  had  a  de- 
stabilizing effect.   Significantly,  this  unstable  result 
could  be  avoided  if  the  pilot  used  both  ailerons  and 
rudder  for  control.   Just  why  the  authors  feel  this 
point  to  be  significant  will  be  discussed  later  in  this 
report. 

Time  histories  of  the  various  angles  and  rates  as 
recorded  on  a  typical  flight  are  included  herein.   In- 
asmuch as  the  mode  of  interest  was  the  Dutch  roll,  the 
authors  felt  the  two  most  significant  time  histories  to 
be  those  of  the  sideslip  angle  ((3)  and  the  roll  rate 
(0).   Accordingly,  the  time  histories  of  these  two 
items  were  studied  in  detail  and  will  be  the  items 
discussed  herein. 
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It  can  "be  seen  from  the  time  histories  of  the  air- 
plane response  to  a  step  aileron  input  (Fig.  11),  that 
the  roll  rate  did  contain  a  considerable  amount  of  Dutch 
roll,  and  further,  that  the  Dutch  roll  motion  was  lightly 
damped.   Also,  it  can  be  seen  that  the  proverse  sense  of 
the  aileron  yaw  resulted  in  the  initial  sideslip  generated 
by  an  aileron  input  being  of  opposite  sign  to  the  roll 
rate  response.   This,  v/as  as  might  be  expected  for  an  air- 
plane with  appreciable  proverse  aileron  yaw  and  only  a 
slight  amount  of  yaw  due  to  roll.   However;  this  initial 
sideslip  response  to  an  aileron  input  v/as  the  opposite 
of  what  was  observed  to  happen  on  the  analog  computer 
simulator.   With  the  benefit  of  hindsight,  it  can  be 
stated  that  this  difference  was  probably  due  to  the  fact 
that  the  effect  of  the  trim  angle  of  attack  (Cyo)  had 
not  been  considered. 

Prom  the  time  histories  of  the  attempts  to  fly  the 
resultant  configuration,  (Pig.  12)  it  can  be  seen  that 
attempts  by  the  pilot  to  control  the  oscillation  resulted 
in  excursions  of  increasing  magnitude  for  both  the  side- 
slip angle  and  the  roll  rate.   Releasing  the  controls, 
(Fig.  13)  on  the  other  hand,  resulted  in  decreasing 
amplitudes  of  motion.   However,  once  large  amplitudes  of 
motion  had  been  obtained,  what  appeared  to  be  a  limit 
cycling  of  the  autopilot  prevented  a  complete  damping  of 
the  open  loop  motion.   An  important  piece  of  information 
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that  can  be  obtained  from  these  time  histories  is  the 
phase  relationship  of  the  sideslip  angle  (  fb   )  and  the 
roll  angle  (  <fi   •)  -  the  latter  being  assumed  to  lag  the 
roll  rate  by  a  quarter  of  a  cycle.   In  particular,  it 
can  be  seen  that  the  sideslip  angle  and  the  roll  angle 
are  almost  exactly  in  phase.   This  was  also  in  contrast 
to  the  phasing  obtained  by  the  electronic  analog  simulator. 
Thus  it  was  noted  that  both  the  analog  simulator  and  the 
NAvion  simulator  exhibited  a  pilot  induced  oscillation 
in  spite  of  the  fact  that  the  character  of  the  motion  was 
virtually  opposite  for  the  two  cases.   That  the  sideslip 
resulting  from  a  step  aileron  input  was  also  different 
for  the  two  systems  suggested  the  possibility  that,  in 
spite  of  the  apparent  difference,  the  same  destabilizing 
mechanism  was  involved  in  each  case.   Accordingly,  the 
aileron  yaw  was  reversed  from  the  proverse  to  the  adverse 
sense  and  the  tests  repeated.   Time  histories  (Fig.  14) 
show  that  although  the  response  of  sideslip  angle  to  an 
aileron  input  was  reversed,  the  percent  of  Dutch  roll 
in  roll  was  maintained  at  a  high  level  and  the  Dutch 
roll  remained  lightly  damped.  However,  this  configuration 
could  be  controlled  by  the  pilot.   This  is  not  to  say 
that  the  resultant  configuration  was  a  good  -  or  even 
acceptable  -  airplane  from  a  handling  qualities  stand- 
point.  The  damping  of  the  Dutch  roll  mode  was  so  slight 
that  unpleasant  oscillations  developed  rather  easily  - 
particularly  if  rapid  control  reversals  were  attempted. 
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However,  it  is  significant  that  the  pilot  was  able  to 
control  the  aircraft,  and,  in  fact,  damp  the  oscillations 
when  they  occurred  using  only  normal  aileron  control 
movements.   This  seemed  to  confirm  the  hypothesis  that 
the  destabilizing  mechanism  for  the  X-15  pilot-aircraft 
combination  (as  manifested  by  the  analog  simulator)  was 
of  a  similar  character  to  what  was  obtained  with  the 
T\ Avion.   In  particular,  it  could  be  seen  that  the  in- 
stability resulted  when  the  pilot,  in  an  attempt  to 
control  roll  angle,  succeeded  in  reducing  the  roll  angle 
only  at  the  expense  of  increasing  the  sideslip  angle. 

During  this  flight  phase,  flight  of  the  basic 
(closed  loop  unstable)  configuration  using  both  ailerons 
and  the  rudder  to  maintain  lateral -directional  control 
was  tried.   The  results  showed  that  by  using  the  rudder, 
the  sideslip  angle  could  be  controlled  and  the  aircraft 
handled  comparatively  easily.   Rapid  reversals  presented 
a  problem  as  it  was  difficult  to  keep  the  aircraft  in 
balanced  flight  throughout  the  maneuver  and  a  rather 
large  Dutch  roll  oscillation  often  resulted.  However, 
by  using  rudders,  the  oscillations  could  be  damped  within 
a  few  cycles,  (Pig.  15).   This  also  tended  to  confirm  the 
importance  of  the  sideslip  angle  to  the  oscillation. 

Aerodynamic  derivatives  other  than  the  aileron  yaw 
were  varied  and  a  pilot  opinion  of  these  variations 
obtained.   However,  even  though  many  of  these  derivatives 
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were  found  to  have  an  important  effect  upon  the  handling 
qualities,  none  were  found  to  "be  so  directly  involved  in 
the  actual  destabilizing  action  of  the  pilot  as  the 
aileron  yaw,  but  rather  were  indirect  contributions  to 
the  problem.   As  such,  the  effects  of  these  other  stability 
derivatives  were  probably  of  importance  in  establishing 
the  conditions  believed  to  be  necessary  for  a  problem 
of  this  type  to  exist.   These  conditions  are  further 
discussed  in  the  following  section. 

The  reason  for  the  difference  in  phasing  of  roll 
and  sideslip  between  the  computer  and  NAvion  is  believed 
to  lie  in  the  inability  to  directly  simulate  the  side 
force  equation  with  the  NAvion.   An  attempt  was  made  to 
adjust  the  NAvion  stability  derivatives  in  such  a  way  as 
to  allow  for  the  effect  of  trim  angle  of  attack  ( <Xo )  • 
Plight  test  of  this  altered  configuration,  however, 
indicated  that  the  phasing  had  not  been  changed. 

The  fact  that  the  NAvion  did  not  yield  the  same 
phasing  of  sideslip  and  roll  as  was  obtained  with  the 
analog  computer  was  not  considered  by  the  authors  to  be 
a  shortcoming  of  the  flight  simulation.   It  was  felt, 
rather,  that  this  difference  in  phasing  helped  to 
emphasize  the  more  general  nature  of  the  physical 
explanation  of  the  instability  mechanism  proposed  here- 
in -  that  is,  since  the  same  character  of  pilot  induced 
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instability  was  manifested  in  both  of  the  two  different 
types  of  systems  investigated,  the  problem  is  not  limited 
to  one  particular  aircraft.   Furthermore,  since  the 
analysis  of  the  mechanism  of  instability  is  consistent 
with  both  of  the  two  types  of  motion  phasing,  the  appli- 
cations of  the  analysis  also  become  less  limited. 

Summary  of  Results; 

Quite  obviously  the  different  phases  of  this  program 
have  resulted  in  a  multi-sided  approach  to  the  basic 
problem  of  a  pilot-induced  Dutch  roll  instability.   In 
general,  it  can  be  stated  that  each  phase  of  the  investi- 
gation contributed  something  unique.   Thus,  in  this  section 
of  the  report,  these  separate  and  distinct  contributions 
will  be  summarized  in  an  effort  to  show  how  each  con- 
tributed toward  the  formulation  of  what  is  believed  to 
be  a  correct  physical  explanation  of  the  destabilizing 
phenomenon. 

Prom  the  root  locus  investigations  discussed  pre- 
viously, it  was  determined  that,  with  a  positive  dihedral 
effect t   an  increase  in  the  trim  angle  of  attack  (  Q£0) 
had  a  destabilizing  effect  upon  the  characteristic  open 
loop  response  of  the  X-15.   Further,  a  study  of  the 
equations  of  motion  showed  that  the  effect  of  increasing 
the  trim  angle  of  attack  was  to  increase  the  coupling 
of  the  sideslipping  and  rolling  motions  of  the  aircraft o 
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In  addition,  a  study  of  the  co-factor  of  the  trim  angle 
of  attack  ((X0)  terra,  in  conjunction  with  the  root  locus 
shown  in  Figure  8c,  indicated  that  the  phasing  of  the 
motions  was  of  such  a  nature  that  destabilizing  moments 
were  generated.   Even  more  signif icantly,  the  root  locun 
analysis  showed  that,  under  certain  conditions,  the  normal 
inputs  of  a  pilot  attempting  to  keep  the  wings  level 
had  a  destabilizing  effect  upon  the  Dutch  roll  mode. 
It  was  then  reasoned  that  inasmuch  as  pilot  response 
time  had  not  been  considered  in  the  analysis,  the  de- 
stabilizing effect  of  the  human  controller  could  not 
have  been  the  result  of  erroneously  applied  corrections 
to  unwanted  deviations  in  roll  angle,  but  rather  had  to 
be  due  to  the  indirect  effects  of  the  resultant  aileron 
deflections.   This  indirect  effect  was  seen  to  be  the 
conversion  of  angle  of  attack  into  sideslip.   Thus  it 
was  concluded  that  the  attempt  of  a  pilot  to  control 
roll  angle  resulted  in  the  generation  of  a  sideslip 
angle.   Further,  the  phasing  between  the  rolling  and 
sideslipping  motions  was  of  such  a  nature  that  one 
motion  tended  to  reinforce  the  other. 

The  results  obtained  as  a  consequence  of  the  analog 
computer-simulator  phase  of  this  program  complemented 
the  theoretical  investigation  and,  in  addition,  provided 
an  explanation  of  the  actual  destabilizing  mechanism. 
The  computer  studies  confirmed  the  importance  of  the  trim 
angle  of  attack  to  the  oscillatory  behavior  of  the  system. 
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Subsequently  the  root  locus  investigation  was  under- 
taken and  the  results  discussed  above  obtained. 

Study  of  the  time  histories  obtained  during  the 
computer  phase  of  the  investigation  showed  that  the 
Dutch  roll  oscillation,  for  the  equations  represent- 
ative of  the  X-15,  was  characterized  by  angles  of  roll 
and  sideslip  which  v/ere  of  opposite  sign,   Further, 
the  open  loop  response  to  a  step  aileron  input  showed 
that  an  input  producing  a  positive  roll  resulted  in  an 
initial  positive  sideslip  (Pig.  9).   This  was  not  what 
was  expected,  inasmuch  as  the  aileron  yaw  was  proverse 
and  the  yaw  due  to  roll  was  small.   However,  as  has  been 
previously  explained,  the  conversion  of  angle  of  attack 
into  sideslip  as  the  result  of  a  rolling  motion  could 
produce  an  effect  of  this  type. 

By  combining  the  two  effects  noted  above  -  which 
would  be  the  natural  result  of  a  pilot  attempting  to 
control  a  Dutch  roll  motion  with  ailerons  -  it  could  be 
seen  that  an  attempt  to  control  the  angle  of  bank  in  a 
normal  manner  would  result  in  an  input  tending  to  in- 
crease the  angle  of  sideslip.   With  the  high  degree  of 
coupling  present  between  sideslip  and  roll,  as  evidenced 
by  the  considerable  Dutch  roll  contribution  to  the  roll 
rate  time  history,  an  unstable  behavior  resulted.   This 
hypothesis  was  checked  by  artificially  increasing  the 
proverse  aileron  yaw  to  such  an  extent  that  the  initial 
sideslip  generated  by  a  step  aileron  input  was  changed. 
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That  is,  an  input  which  resulted  in  a  positive  roll 
yielded  a  negative  sideslip.   Figures  16  and  17  show 
the  time  histories  resultant  from  this  configuration. 
With  the  increased  proverse  aileron  yaw,  the  configuration 
could  be  controlled  in  spite  of  the  fact  that  the  Dutch 
roll  mode  remained  lightly  damped  and  the  percentage  of 
Dutch  roll  in  roll  remained  appreciable.   The  lightly 
damped  Dutch  roll  mode  did  result  in  easily  induced 
oscillations  but,  unlike  the  basic  configuration,  pilot 
attempts  to  control  roll  angle  did  not  destabilize  the 
system.   In  fact,  as  might  be  expected,  it  was  found  that 
attempts  to  control  the  roll  angle  in  a  normal  manner 
also  helped  to  reduce  the  excursions  of  the  sideslip 
angle. 

A  by-product  of  the  computer  study  was  the  confirm- 
ation that  the  frequency-  of  the  oscillation  was  not  a 
significant  contributor  to  the  basic  instability.   Quite 
obviously,  this  does  not  mean  that  frequency  can  always 
be  dismissed  as  insignificant,  but  only  that  an  instability 
of  the  type  studied  could  exist  even  at  frequencies 
completely  within  the  control  limits  of  a  human  being. 
In  this  phase,  it  was  found  that  motions  with  a  period 
as  high  as  five  seconds  were  completely  unmanageable 
as  long  as  a  normal  pilot  technique  v/as  used. 

The  flight  phase  of  this  investigation  yielded 
results  which  essentially  complemented  those  obtained 
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from  the  fixed  base  simulator.   This  was  in  spite  of  the 
fact  that  the  character  of  the  oscillatory  motion  obtained 
with  the  NAvion  variable  stability  airplane  was  consider- 
ably different  from  what  was  observed  with  the  fixed 
base  simulator. 

As  outlined  previously,  the  connection  seemed  to  be 
that,  in  either  case,  an  attempt  by  a  pilot  to  control 
the  roll  angle  by  manipulation  of  the  ailerons  resulted 
in  an  unfavorable  generation  of  sideslip.   Inasmuch  as 
the  rolling  and  sideslipping  motions  were  highly  coupled, 
the  increased  sideslip  generated  additional  rolling 
moments  which  in  turn  required  even  greater  control 
movements.   The  result  was  a  coupled^  rolling,  sideslipping 
motion  of  ever  increasing  amplitude.   That  the  oscillations 
could  be  contrdlled  -  or  even  damped  -  by  judicious  use  of 
the  rudder,  substantiated  the  validity  of  the  foregoing 
explanation  of  the  instability  mechanism. 

The  result  of  the  combined  studies  was  that  certain 
combinations  of  roll  angle,  sideslip  angle,  and  side- 
slip resultant  from  an  aileron  deflection,  were  found 
to  be  unfavorable.   In  particular,  it  was  noted  that 
the  pilot-induced  unstable  behavior  occurred  only  when 
the  pilot,  in  an  effort  to  return  to  a  wings  level 
condition,  could  not  reduce  the  bank  angle  without 
increasing  the  sideslip  angle.   Thus  an  "effective" 
proverse  aileron  yaw  was  found  to  have  a  destabilizing 
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effect  if  the  Dutch  roll  oscillation  was  of  such  a  nature 
that  the  bank  angle  and  the  sideslip  angle  were  in  phase. 
For  an  "effective"  adverse  yaw,  the  opposite  phasing  was 
unfavorable.   Multiple  exposure  photographs  depicting 
the  two  motions  described  above  are  included  herein. 
(Pigs.  18  and  19).   The  term  "effective"  as  used  refers 
to  the  sense  of  the  initial  sideslip  produced  by  an 
aileron  deflection.   This  initial  sideslip  depends  on 
the  combined  effects  of  Cn&a  r^np  and  O£0f   rather  than  the 
actual  sign  of  the  yawing  moment  produced.   Thus  it  might 
be  termed  an  integrated  effect. 

In  addition  to  the  above  stated  phase  relationships, 
certain  conditions  were  necessary  before  the  instability 
could  occur.   First,  it  was  necessary  that  the  lateral- 
directional  motions  be  highly  coupled.   Physically  this 
was  indicated  by  the  large  oscillations  which  occurred  in 
the  time  histories  of  roll  rate  response  to  a  step  aileron 
input.   Secondly,  it  was  necessary  that. the  Dutch  roll 
mode  be  lightly  damped.   No  formal  proof  of  this  last 
statement  has  been  offered  herein,  but  it  seems  quite 
obvious  that  a  highly  damped  mode  would  virtually  never 
demand  pilot  attention.   Practically  speaking,  the  pilot 
would  probably  never  be  aware  of  its  existence.   Accord- 
ingly, one  could  hardly  expect  the  pilot  to  ever  initiate 
the  destabilizing  mechanism  described  above.  From  a  root 
locus  standpoint,  an  increase  in  the  Dutch  roll  damping 
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Dutch    Roll     Mo+ion 
ft  and     <p    of     Opposite     Phase 


Fig.     19 
Dutch      Roll     Motion 
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would  shift  the  Butch  roll  poles  and  zeroes  (Pigs.  8a 
and  8b)  to  the  left.   This  would  obviously  have  a 
stabilizing  effect. 

From  the  discussion  of  the  destabilizing  mechanism, 
it  is  apparent  that  the  phasing  between  the  bank  angle 
( <p  )  and  the  sideslip  angle  ( ft )  is  an  important  indicat- 
or of  how  a  closed  loop  instability  might  develop.   In 
particular,  this  phase  relationship  determines  the  sense 
of  the  effective  aileron  yaw  that  can  be  regarded  as 
unfavorable.   Accordingly,  a  study  was  undertaken  in  an 
effort  to  ascertain  what  parameters  determined  this  phase 
relationship.   This  study  is  outlined  in  Appendix  C. 
Briefly,  the  phase  angle  between  <p   and  fd  was  found  to 
be  heavily  dependent  upon  the  sign  of  the  termr-^Oflft-y-S 
This  same  combination  of  derivatives  has  been  found  to  be 
important  in  the  determination  of  the  ^^  parameter0 

Further  analysis  led  to  the  conclusion  that  the 
phasing  of  the  Dutch  roll  in  roll  could  be  expected  to 
bear  some  relationship  to  the  character  of  the  roll  rate 
response.   Figures  9  and  16  showed  that  a  distinct 
difference  existed  between  the  characteristic  time 
constants  of  the  roll  mode  for  the  favorable  and  un- 
favorable situations.   In  particular,  the  time  difference 
was  seen  to  be  approximately  equal  to  one  half  of  the 
Dutch  roll  period.   This  effect  is  discussed  in  detail 
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in  Appendix  D.   Briefly,  a  shift  of  180  degrees  (Dutch 
roll  in  roll)  was  seen  to  coincide  with  a  shift  through 
unity  of  the • magnitude  of  the  J&T  parameter. 

In  summary,  for  the  type  of  pilot  induced  instability 
outlined  herein,  the  arguments  presented  appear  to  "be  a 
physical  way  of  describing  the  phenomenon  depicted  by 
the  root  locus  studies  of  Taylor.   The  theoretical 
effort  conducted  in  conjunction  with  this  investigation 
seems  to  show  that  the  onset  of  an  unfavorable  combination 

of  phase  relationships  is  coincident  with  the  parameter 

CO* 

described  as  -j§J  becoming  greater  than  unity. 

Although  only  a  limited  number  of  cases  of  pilot- 
induced  instabilities  were  studied  in  detail  in  con- 
junction with  this  investigation,  the  authors  have  seen 
no  reason  to  suspect  that  the  destabilizing  mechanism 
as  outlined  herein  cannot  be  regarded  in  a  broader  sense. 
Thus  it  is  the  opinion  of  the  authors  that  an  explanation 
of  this  type  might  well  be  valid  for  other  instabilities 
known  to  be  a  function  of  pilot  gain. 


-68- 
SBCTIOK  V 

COKCLUSIOKS 

As  a  result  of  this  investigation,  the  following 
conclusions  were  drawn: 

1.  The  NAvion  variable  stability  aircraft  can  be 
used  for  a  general  study  of  pilot-induced 
lateral -directional  instabilities. 

2.  A  pilot  induced  instability  of  the  type  investi- 
gated can  be  successfully  demonstrated ,and 
analyzed  with  a  fixed  base  simulator  incorpo- 
rating an  electronic  analog  computer. 

3.  For  this  type  of  pilot  induced  instability  of 
the  Butch  roll  mode  to  occur,  certain  conditions 
were  found  to  be  required.   These  are  listed  in 
a  qualitative  fashion  as: 

a.  A  lightly  damped  Dutch  roll  transient; 

b.  Highly  coupled  rolling  and  sideslipping 
motions  -  as  manifested  by  a  large  per- 
centage of  Dutch  roll  in  roll. 

4.  If  the  conditions  listed  in  item  (3)  were  prevalent, 
certain  phase  relationships  between  the  roll  angle, 
sideslip  angle ,  and  the  effective  aileron  yaw  were 
found  to  result  in  an  unstable  closed-loop  system. 
In  particular,  if  the  Dutch  roll  oscillation  showed 
the  roll  angle  and  the  sideslip  angle  to  be  of  the 
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same  sign,  an  unstable  system  resulted  if  the 
effective  aileron  yaw  was  proverse.   For  the 
opposite  phasing,  an  effective  adverse  aileron 
yaw  resulted  in  an  unstable  closed  loop  system. 

5.  For  the  type  of  pilot  induced  instability 
studied,  the  frequency  of  the  oscillation  was 
found  to  be  an  unimportant  parameter. 

6.  Even  a  marginally  stable  system  (open  loop) 
could  be  managed  quite  easily  if  a  direct  means 
for  controlling  the  sideslip  angle  were  provided. 
Accordingly,  a  system  which  was  closed  loop  un- 
stable with  only  aileron  control  was  found  to  be 
manageable  if  the  rudder  was  usedo   Also,  the  & 
technique  could  be  used,  although  greater  care 
was  required. 

7.  It  was  further  concluded  that  this  study  provided 
a  physical  interpretation  of  the  type  of  pilot- 
induced  instability  predictable  by  the  method 
described  by  Taylor  (Ref.  2). 
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TABLE  I 
PHYSICAL  CHARACTERISTICS  OF  THE  NAvion 


WING  DATA 

Total  Area  (including  ailerons,  flaps 

and  19.87  ft.  covered  by  fuselage) 

Span 

Aspect  Ratio 

Taper  Ratio 

Dihedral  Angle 

Root  Chord 

Mean  Aerodynamic  Chord 

Incidence  Angle 

Root 

Tip 
Sweepback  of  Leading  Edge 
Twist 

Geometric 
Airfoil  Section 

Root 

Tip 
Flaps,  40  deg. ,  Plain 

AILERON  DATA  (For  one  aileron) 
Area 
Span 


183.34  ft. 
33.38  ft. 
6.05 
0.54 

7.50  deg. 
7.20  ft. 
62.35  in. 

2.0  deg. 

-1.0  deg. 

3.0  deg. 

3.0  deg. 

KACA  4415R 
NACA  6410R 


2.16  ft.' 
61.99  in. 


Deflection 


30  deg.  up;  20  deg.  dwn. 
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TABLE  I  (continued) 

AILERON  DATA  (continued) 
Control 

Aerodynamic  Balance 
Static  Balance  (Outboard  end  of 
each  aileron) 
Trim  Tab  (Right  aileron) 
Ratio  of  Aileron  Chord  to 
Wing  Chord 

HORIZONTAL  TAIL  DATA 

Total  Area  (Including  2.37  ft.2 

covered  by  fuselage) 

Span 

Aspect  Ratio 

MAC 

Airfoil  Sections,  Root  and  Tip 

Incidence  Angle 


ELEVATOR  DATA 
Total  Area 
Span 

Deflection 
Deflection  Trim  Tabs 
(32  in.  span,  4-1/2  in.  Chord) 
Root  Chord 
Tip  Chord 


Wheel  Throw 
Prise-type  nose 

Streamlined  weight 
Fixed  Bend  Tab 

0.284 


43.05  ft.' 
13.17  ft. 

4.02 

3.34  ft. 
NACA  0012 
-3.0  deg. 


14.10  ft. 
73.58  in. 
30  deg.  up;  20  deg.  dwn. 


+30  deg. 
1.5  ft. 
1.0  ft. 
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TABLE  I  (continued) 

VERTICAL  TAIL  DATA 

Area  (Including  2.57  ft.2  blanketed  by 

fuselage  and  excluding  1.84  ft.   dorsal 

fin  area)  12.93  ft.2 

Span  4.05  ft. 

Airfoil  Section 

Root  NACA  0013.2  Modified 

Tip  NACA  0012-64  Modified 

Incidence  Angle  (With  respect 

to  FRL)  2  deg.  Nose  Left 

RUDDER  DATA 

Area  6.05  ft.2 

Deflection  17  deg.  L;  23  deg.  R. 
Trim  Tab  Fixed  Bend  Tab 

Rig  (Angle  with  respect  to 
fin  center  line)  3  deg.  right 

FUSELAGE  DATA 

Fuselage  Length  Over-all  27.25  ft. 

Width,  Maximum  4.14  ft. 

Depth,  Maximum  4.40  ft. 

Fineness  Ratio  6.2 

POWER  PLA1\T  DATA 

Airplane  is  powered  by  one  Continental  E-185  engine. 
Maximum  continuous  rated  horsepower  at  sea  level 
185  at  2300  RPM. 


100 
86  in. 
26.5  deg. 
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TABLE  I  (continued) 

POWER  PLANT  DATA  (continued) 

The  propeller  is  a  Hartzel  hydro-selective  propeller 
with  the  following  characteristics: 

Activity  Factor 

Diameter 

Pitch  (.75R) 

MISCELLANEOUS  DATA 
Weight 
Basic 

Fuel  (40  gal.) 
Pilots  (2  with  parachutes) 
Gross  Weight 
Center  of  Gravity  Position 
Tail  Length 

Autopilot  Mechanical  Gear  Ratios 
Deg. Elevator  per  Deg.  Servo-drum 
Deg. Aileron  per  Deg.  Servo-drum 
Deg.  Rudder  per  Deg.  Servo-drum 


2129 

lbs. 

240 

lbs. 

410 

lbs. 

2779 

lbs. 

29, 

.5$  MAC 

16, 

,88  ft. 

0, 

.919 

0, 

,910 

0, 

,330 
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TABLE  II 
SPECIFICATIONS  OF  INSTRUMENTATION  COMPONENTS 


RATE  GYRQS 
Manufacturer 
Model  Number 
Fower  Input 
Power  Consumption 

Rotor  Speed 

Weight 

Maximum  Turn  Rate 

Potentiometer  Resistance 
Potentiometer  Excitation 


SIDESLIP  VANE 
Manufacturer 
Model  Number 
Potentiometer  Resistance 
Potentiometer  Active  Angle 

Potentiometer  Excitation 


Minneapolis-Honeywell 

JG  7005A-24 

115  volts, 400  cycle  A.C. 

32  watts  (starting), 

13  watts  (running) 

20r000  RPM 

1.75  lbs. 

Modified  to  approximately 

30°/sec. 

530  ohms 

Gyro  #9  -  24  volts  A.C. 

Gyros  #10,  #11  and  #13  - 

30  volts  A.C. 

Giannini 

2516 

1997  ohms 

Yaw  channel  -  +  30° 

Roll  channel-  +45° 

Yaw  channel  -  15  volts  A.C. 

Roll  channel-  30  volts  A.C. 
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TABLE  II  (continued) 

TELEMETER  TRANSMITTER  UNIT 
Manufacturer 
Model  Number 
Input  Signals 
Information  Channels 


Sampling  Rate 

Accuracy 

RP  Power  Output 

Frequency  Range 

Frequency  Stability 

Primary  Power  Requirement 

Weight 


AS  GOP 
DT-4 

0-5  volts  D.C. 
43  plus  two  for  synchro- 
nization of  ground  station 
20  RPS 

+  ifo   of  full  scale 
4  watts 
215-235  MC 
+  0.05$ 
28  volts  D.C. 
19.6  lbs. 
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TABLE  III 

Dimensional  derivatives  of  the  X-15  at  Mach  =  3.0,  cxo»  10° 

Lp   =  +9.02  Np   =  -.00451 

Lp   =  -.212  Nr   =  -.0451 

Lr   =  +.172  N^q   =  +.454 

Np   =  +2.75  L^   =  +8.40 

TABLE  IV 
Kon  dimensional  stability  derivatives  of  the  X-15 
at  Mach  =  3.0,  (X0 ■  10° 

%  -  +.044 

Cn  =  -.280 

Co  =  +.228 

Cn^  =  +.315 

Cnp   =   -.140 

TABLE  V 

Autopilot  potentiometer  feedback  gain  constants. 

k 

4  =  +.541  deg.  rudder/deg.  sideslip 

5  =  -.121  deg.  rudder/deg. /sec  roll  rate 
kg     =  -.111  deg.  rudder/deg. /sec  yaw  rate 

kg     =  +.609  deg.  aileron/deg.  sideslip 

kg     =  -.0756  deg.  aileron/deg. /sec  yaw  rate 

k7^k13  =  -.345  deg.  rudder/deg.  aileron 


CnP       = 

-1.40 

Cn&a     = 

+  .052 

c*6a    ■ 

+  .041 

Cu       = 

+  .440 

%     '= 

-1.38 
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A3J-ENDIX  A 
Gale ulat ions  for  the  Plight  Simulation: 

Non-dimensional  stability  derivatives  chosen  as 
typical  for  the  X-15  at  M  =  3.0,  ex  =  10°  are: 


=  -1.38 


Cj^    .  +   .044 

Cnp    =   -    .140 

cH    = 

Cflp    =   -    .280 

Cnr   =  -  1.40 

c^ 

Cj?p    =   +    .228 

C^a=   4    .052 

Cn/3    =  +    .315 

Cfl6a  =  +    .041 

Por  one  -  g 

flight: 

« 

f- v-: 

1/2   CL  V^   S 

13,445 

1/2    (.44X3  x  1005)2(200)    =   3.36  x  10~5lb-sec2/ft4  (1) 

Where  1005  ft/sec  is  the   speed   of   sound  corresponding  to 
an  altitude  for  which  p   =  3.36  x  10~5lb-sec2/ft4  (2) 

q  =  1/2  pV2   =  1/2   (3.36xlO~5X9.1xl06)   =  153  lh/ft2  (3) 

T  =  M/PSV  =  3MiT0-5   (200)(3015X32.174)=  20'6  (4) 

/*=  M/i°Sb   =  TTfelO"5   (200)(32. 174X22. 36)    =  2790  (5) 


2  ~      ~        2 

=   3,348  - 

137W 


k  *   =  Ix  =  3,348  x  32.174   =  8.02  ft^  (6) 


k^2   =   I-   =   78,691  x  32.174   =  188.4  ft2  (7) 

z         „£       T51445 

J      =   2fkx^|2   =  2      8.02   =    .0321 

x        Vb~7  5^T 


(8) 


J     =  2    /kz]2  =  2     188.4  =   .753  (9) 
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Values   for  dimensional  derivatives  were   obtained  as 
follows: 

Lft  =  i^A=    .044(2790)        =   +   9.02  (10) 

Lp  =£M-  =  --2Q0  ■  --212  (ID 

P      2TJ*         2(20. 6)(. 0321) 

Lr  =  S^T  =    _    -228  =    +.172  (12) 

2^J*       2(20. 6)(.  0321) 

Nft  =/ATr^=  2790(.31g)     =  +  2.75  (13) 


eTJ£     2  <26.6X.724> 


Np  =;^^-=  "   *1* =  -.00451  (14) 


2-TJt     2   (20.6)(.?24) 

i 

K^^^r2  ■    .052(2790)        =   +.454  (16) 
a   T'Jt      424.4  (.753) 

Ua=  ^2^=    .041    (2790)        =   ♦   8.40  (17) 
T\J*     424.4  (.6321) 

The   equations  which  express   the   relationship  "between  the 
simulator  airplane  derivatives   and   the   simulated  deriva- 
tives are  derived  in  Reference  4.      The   resultant  express- 
ions for  the  autopilot  feedback  gain  constants  are: 
k4   =  1.312   -  2.415  x  10"3Lj3     -  0.2774  Nip  (18) 
k5   =   -0.123  -  2.415  x  10"5_p    -  0.2775  Np  (19) 
k6   =   -0.124   -  2.415  x  10~3Lr    -  0.2775  Nr  (20) 
kQ   =  0.346   +  0.0490  L^    +  0.0278  Nfc  (21) 
kg   =  0.373   +  0.0490Lp+  0.0278  Np  (22) 
k1Q   =  0.0828  +  0.0490  Lr   +  0.0278  Nr  (23) 
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k7  =  -  20.5(N6g/L<W)  +  0.178  (24) 

k13     2.05CN^Q/L6a)^3.^^ 

Values  obtained,  from  the  foregoing  equations,  for  feedback 
gain  constants  were: 

k.  =  +  .541  deg.  rudder/deg.  sideslip 

kr  =  -  .121  deg.  rudder/deg/aec  roll  rate 

kg  =  -  .111  deg.  rudder/deg/sec  yaw  rate 

kg  =  +  .609  deg.  aileron/deg  sideslip 

kq  =  +  .362  deg.  aileron/deg/sec  roll  rate 

k10  as  -  .0756  deg.  aileron/deg/sec  yaw  rate 

ky/k^,  =  -  .345  deg.  rudder/deg.  aileron 

Cockpit  potentiometer  settings  were  obtained  from  the 
calibration  curves  for  these  values  of  feedback  gain. 
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APPEUDIX  B 


The  following  is  a  derivation  of  the  root  locus 
equation  used  to  study  the  effect  of  pilot  gain. 

A  pilot  response r  proportional  to  roll  angle,  of 
the  form  K(l  +  .57d)  is  assumed. 

The  determinant  setf  in  non-dimensional  form,  Is 
as  follows: 


t 

d+ 

* 

6a 

C^-Zd 

-z 

CL+  2c<0d 

0 

/aQJj* 

QVz. 

cj£d~J*dz 

MC^a 

/*£"/& 

QVz~ 

J£d 

Q^Bd 

MC"6q 

0 

0 

K(l    +    -5Td) 

-I 

(1) 


Expansion  of  the  above  can  be   symbolized   in  determinental 
form  as: 


-Z 

Cflr/2. 


Cl  +    ZotQd 


+ 


K(l  +  .57d) 


CV/5-Zd 


-a 

On     /.     - 


■  nr. 


72 


J*d 


O 
uCn<3Q 


-  0 


(2) 
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Expansion  of  the  fir^t  determinant  set  will  yield 
a  quartic  equation  whioh  can  be  divided  "by  2J  J  and  then 
factored  into  a  form  from  whioh  the  characteristic  modes  of 
transient  motion  can  be  deduced.   Thus  this  is  the  character- 
istic equation  and  will  be  symbolized  as: 

|D|  '-   (X  +  XR)(X  +  Xs)(Xa  +  2.^C0ny+  OOn^2")         0) 

The  remaining  determinant  can  be  expanded  and  the 
sum  represented  as  follows: 

1  d|  -  K(i+.57d)  Lcfejcy  &)(  %  -  Jid)  +  ^/xCn^ 

-yuCniaLfcvp-^)^  ^mC^JJ 
This  can  be  rearranged  to  yield: 


(5) 


|D|+  K(l  +  57d)/ACjt&,  j[dfcl  fcnr+  J^CYfi)dl  +°^L^nel 

>Jx.     (  2JX  2.  J  * 

dux,  zjs  L     v  -&    ^ciTVJj 

The  terms  multiplying  the  human  transfer  function  are 
seen  to  be  a  quadratic  which  is  known  to  yield  a  complex 
pair  of  roots.  Accordingly  the  above  may  be  symbolized 
as  follows:   (d  is  replaced  by  A  ) 

|D|  +^5Ma.K(l  +  .57A)(XL+  2.^COn0  +  COn/)        (6) 
Jx 

The  root  locus  equation  is  thus  found  to  be: 

(Xz  +  Z%<p'odn<p  X  +  con^2-)(l  +  .57 X)  Jx 
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APl'ENDIX  C 


Discussion  of  Phasing  of  goll  and  Sideslip; 

One  of  the  parameters  discussed  in  conjunction  with 
the  instability  studied  herein  is  the  phase  relation  of 
roll  angle  and  sideslip  angle.   This  relationship  is  the 
subject  of  the  following  discussion. 

Prom  the  basic  determinant  for  the  lateral- 
directional  motion  of  the  airplane,  given  here  in  Laplace 
operator  form  as: 


/^Cn/o 


-z 

Z 


Cn. 


-  J*S 


7L 


6a 
0 

_  LxCnSa 


(l) 


one  may  form  the  ratio  <p  to/5  as: 


_£ 
P 


O 

_  uCn6a 
1   '  s 


-z 

2. 


-  JxS 


-z 

2. 


O 


7  s 


Cn£5 


(2) 
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This  may  be  written  as: 


Vp    "  a2  *  A  s  +  B  (5) 


C   8      +   D  s   +  E 


where 

A         Cn6a    Cjlr 
CJlSq     ZJ£ 

-^ 

Cnr 
9    1 

B   =    i^L  +    St 

fiCnr   _ 

Cn^aCfir 

G  _     2.JX  Cn6a    _    2.o(o  J  a. 
Ci^a 

D  =  CnP  -  ClJ£  +  a0Cnr  -  fo&q  (CCp  4-  ocoCjid 
*    ClCht   __   Cn^aGuCfir 


For  the  condition  of  low  dynamic  pressure  (153  lb/ft  ), 
with  cyG  =  10°  this  ratio  can  be  written 

—  s     »2  +  1.81s  +  962.1         rt„  ,.. 
5 or           (.4; 

P  -,182s*  -  ,366s  -  .3716 


0   = 


b  +  _.9i_±3i.oj,     Constant 


s  +  1.01  +  .957  j 

These  roots  are  plotted  on  the  complex  plane  and  yield  a 
phase  relation  for  0/3  of  -173  degrees  (Pig.  20). 

It  appears  from  an  examination  of  the  above  expression 
of  ^ / fb    that  in  those  cases  where  the  Dutch  roll  mode  is 
lightly  damped,  the  damping  terms  of  the  numerator  and 
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D 


FIG.  20 

<t>/p    PHASE 
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real 
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<t>//3    PHASE   =    Zh  *  ZB  -  /c  -  /o 

=    -173° 
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denorainator  will  also  be  small.   This  means  that  all  the 
roots  will  lie  close  to  the  imaginary  axis.   Further, 
by  analyzing  the  plot  of  these  roots,  it  is  evident 
that  the  phasing  of  the  0/A  motion  is  primarily  de- 
pendent upon  the  relative  position  of  the  numerator 
root  and  the  Butch  roll  root.   If  the  numerator  root 
has  an  imaginary  part  greater  than  the  imaginary  part 
of  the  dutch  roll  root,  the  phasing  between  0  and  jd 
will  be  approximately  180  degrees.   If  the  Butch  roll 
imaginary  part  is  the  larger,  then  the  phasing  will  be 
nearly  zero  degrees. 

Those  parameters  which  affect  the  imaginary  pr:.rt 
of  the  Butch  roll  root  are  difficult  to  isolate  since 
all  approximations  applied  to  the  condition  under  investi- 
gation yielded  poor  results.   However,  it  is  evident 
that  a  major  contribution  comes  from  /M^n/°.   This 
combination  of  parameters  also  appears  in  coefficient  B 
in  the  expression  for  0//3  ,  so  it*  influence  in  shift- 
ing the  two  roots  is  equal.   The  only  other  significant 
term  in  coefficient  B  is  -Z.uGOft  ^n6q   .   it  is  through 
this  expression  that  a  major  shift  in  the  numerator  root 
is  achieved.   In  the  case  under  consideration,  this  one 
term  is  large  enough  to  shift  the  numerator  root  below 
the  Butch  roll  root,  changing  the  phasing  to  zero  degrees. 
This  shift  would  occur  if  one  of  the  three  parameters  had 
a  sign  change. 
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This  is  merely  additional  supporting  evidence  for 

the  argument  that  the  phase  relation  "between  <p   and  fb 

t 
and  the  yaw  due  to  aileron  are  important  parameters 

in  the  explanation  of  pilot  induced  lateral-directional 

instabilities.   It  also  points  out  the  problem  which 

will  arise  in  any  airplane  -  pilot  combination  where 

C|},5  changes  sign  within  the  performance  envelope.   If 

this  change  of  sign  occurs  in  an  area  where  other 

conditions  are  favorable  to  pilot-induced  instabilities, 

then  these  instabilities  may  well  be  encountered. 
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APPSNDIX  D 
Discussion  of  Dutch  Roll  to  Roll  Mode  Phasing: 

In  an  attempt  to  shed  more  light  on  the  mechanism 
producing  the  type  of  instability  under  consideration, 
the  following  analysis  of  Dutch  roll  to  roll  mode  was 
made.   Of  particular  interest  was  the  phase  relation 
of  the  two  motions. 

Again  the  basic  determinant  for  the  lateral- 
directional  motion  of  the  airplane  is  given  here  in 
Laplace  operator  form. 


/ACifo 


-z 

CQr 
2. 

Z 


0 


Forming  an  expression  for  <p   in  terms  of  this 

determinant  and  applying  a  step  input  yields  the  follow- 

ing  determinant  for  (f>  . 

-Z         O 

Or/Z  -fxCjlSa 

Qgr-Jts     -yuCn^q 


cp    a 


GYp-Zs 


This  may  be  written  as: 


-Z 


^5 
Z 


Qnr_J    s  CnpS 

Z  Z 


<P     * 


A( 


Bs 


C 


(s-X.Xs-X^s-X^s-X^ 
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where 

C£6q 


A 


/^ 


R    Qic_  Cn6g     (Cnr  +  Jz-Cvft) 

and  the  denominator  represents  the  characteristic  of 
motion. 

For  the  condition  of  low  dynamic  pressure  (153  lb/ft  ) 
with  CYo  =  10°,  this  expression  becomes: 

0  _  3565  (s^4-  1.814-s  +  9b  P) 
(s-X1Xs-\lXs-X3Xs-X^ 

The  roots  of  the  above  are  plotted  on  the  complex 
plane  (Pig.  21).   It  can  be  seen  that  the  phase  angle 
changes  by  180  degrees  whenever  the  numerator  root  moves 
inside  the  Dutch  roll  root.   Further t   a  comparison  of  the 
roll  traces  (Figs.  9  and  16)  shows  the  same  result. 
Figure  9,  with  effective  adverse  yaw,  shows  a  time  delay 
of  1.5  seconds  to  reach  the  first  peak  on  the  roll  trace. 
Figure  16,  with  the  same  conditions  except  effective 
proverse  yaw,  shows  a  time  delay  of  .5  seconds.   This 
difference  in  time  corresponds  to  one  half  the  period  of 
the  Dutch  roll  motion.   Here  again,  as  in  the  (f>   to  /£> 
relation,  we  observe  a  change  in  phase  of  180  degrees  at 
the  same  point  where  u-V/lOij)  becomes  greater  than  1. 
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